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Abstract 

We study the exact Bethe Ansatz solution of the p + ip Hamiltonian in a form whereby quantum 
numbers of states refer to hole-pairs, rather than particle-pairs used in previous studies. We hnd an 
asymmetry between these approaches. For the attractive system states in the strong pairing regime take 
the form of a quasi-condensate involving two distinct hole-pair creation operators. An analogous feature 
is not observed in the particle-pair picture. 


1 Introduction 

The Px + ipy-wave pairing (or simply p + ip) Hamiltonian rose to prominence through the influential work of 
Read and Green [1] in identifying topological properties of superconducting systems. An exact Bethe Ansatz 
solution of the Hamiltonian appeared in [2], a result which subsequently generated several studies in the 
exactly solvable framework [3-9]. Unlike its counterpart the exactly solvable s-wave pairing Hamiltonian, 
also known as the Richardson model [10,11], the p + ip model exhibits quantum phase transitions. A key 
to gaining a complete understanding of the model’s properties is to understand the distribution of the roots 
of the Bethe Ansatz equations. For example, it is well-established [4, 5] that there is a duality between 
states in the regimes known as weak pairing and strong pairing. There are states in the weak pairing regime 
characterised by the presence of particle-pairs which carry zero energy, and corresponding dual states in 
the strong pairing regime with the same energy. Particular groupings of particle-pairs should be viewed as 
bound, and only arise when certain constraints are satisfied. The existence of these zero energy particle-pairs 
may be interpreted as a form of quasi-condensation. 

Our objective is to revisit the Bethe Ansatz solution of the p + ip model from the hole-pair perspective 
rather than the particle-pair perspective. This involves working with a very closely related, yet distinct, 
second form of Bethe Ansatz solution. There are a couple of motivations for taking this approach. The 
first is that Hirsch has long advocated that “electron-hole asymmetry is the key to superconductivity” [12]. 
As will be established, there is indeed an asymmetry between hole-pairs and particle-pairs which manifests 
in the character of the roots of the associated Bethe Ansatz equations. The second motivation is that 
there have been some very interesting studies of the Bethe Ansatz solutions for pairing Hamiltonians which 
exploit the existence of two forms of solutions. Pogosov and collaborators have determined sets of relations 
between the roots of the Bethe Ansatz equations in the s-wave [13] and “Russian Doll” [14] models, which 
in particular has led to a formula for the ground-state energy in the s-wave case. Also Faribault and 
collaborators [15,16] have used the existence of the two forms of exact solution to facilitate the calculation 
of wavefunction overlaps and scalar products for a general class of systems, which includes the s-wave model. 
Subsequently Claeys et al. [9] have generalised this latter work to accommodate the p + ip model. 

When the hole-pairing picture is adopted, we find that zero energy hole-pairs arise which characterise 
the same duality that was mentioned above. We also find that in addition there exist states characterised 
by infinite energy hole-pairs, in a manner that the sum of the infinite energies is finite. Such infinite energy 
solutions have been previously observed [6,8], but so far a systematic investigation of them has not been 
undertaken. Our main finding indicates a second form of duality, which relates the energies of eigenstates 
of the attractive pairing Hamiltonian to energies of eigenstates of the repulsive pairing Hamiltonian. It also 
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provides a new perspective on the different regimes of the system in terms of hole-pair quasi-condensation, 
which we will discuss later. 

In Sect. 2 we present the Hamiltonian and the known Bethe Ansatz solution in terms of particle-pair 
quantum numbers. We then determine the second solution in terms of hole-pair quantum numbers. Sect. 3 
is devoted to the discussion of dualities. We commence by recalling the duality previously discussed in [4,5]. 
We then continue to establish a second form of duality, and ultimately a third which is a combination of 
the two. Sect. 4 examines these dualities in the framework of the Bethe Ansatz equations. In Sect. 5 we 
undertake numerical solution of the Bethe Ansatz equations, which leads into an investigation of the phase 
diagram in Sect. 6. Sect. 7 briefly compares the results of mean-field approximations with the results of 
the exact solution. Concluding remarks are given in Sect. 8. 


2 The Hamiltonian and exact solution 

We first introduce the Hamiltonian of the pairing model. We take the canonical (i.e. particle number¬ 
preserving) Hamiltonian whose mean-field approximation leads to the Bogoliubov-de Gennes equations 
with order parameter having px + ipy-wave symmetry up to quadratic approximation. Letting Ck, denote 
annihilation and creation operators for two-dimensional fermions of mass m with momentum k = {kx,ky), 
the Hamiltonian reads [2] 

= ikx+iky){k'x-ik'y)clcl^^c-i,>ck,, 

k k^^itk' 

where G is a dimensionless coupling constant which is positive for an attractive interaction and negative for 
a repulsive interaction. For any unpaired fermions the action of the pairing interaction is zero and we can 
thus decouple the Hilbert space into a product of paired and unpaired fermions states, for which the action 
of the Hamiltonian on the space for the unpaired fermions is diagonal in the number operator basis. This 
is known as the Mocking effect [11] and permits an analysis of a simplified version of the Hamiltonian. 

We set Zk = |k| and kx +iky = |k] exp(f())k). It is convenient to introduce the following phase-dependent 
Cooper pair (or hardcore boson) operators b]J = exp(i())k)C]jj.c[Lk = ^-k^, &k = exp(—f(()k)c_kCk = 6-k, and 
set fVk = b]Jbk = A^-k- Using integers to enumerate the unblocked pairs of momentum states, and setting 
m = 1, the Hamiltonian takes the form 

L L L 

HiG) = Y,zlN,-GY,Y. ziZkbjbk. 

k—1 1—1 k^l 

The hardcore boson operators satisfy the following commutation relations 

hA] = - 2N,), [b„bk] = [b], bi] = 0 ( 1 ) 

where / denotes the identity operator, as well as the relations = 0, = Nj. The Hamiltonian may 

be expressed in the compact form 

H{G) = (1 + G)Ho - GQ^Q (2) 


where 


L L L 

ho = y, = E ^ = E 

1^1 1^1 

For later use we note the commutation relation 

L 

[Q\Q] =2Ho-Y,zp 

i=i 


(3) 


(4) 
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which follows from (1). 

The Hamiltonian (2) is our principal object of study. For each solution of the coupled equations 

G-1+2M-T-1 ^ 1 ^ 2 , ^ 

vk ^ 2/fe - zf ~tky^~ 

there is an eigenstate of (2) with energy eigenvalue given by 

M 


E = {i + G)J2yk- 




The eigenstate has the form 


( 5 ) 


( 6 ) 


M 

m = Y[C{y,)\0) 

i=i 

where 

j=i y 

Above, M is a quantum number which denotes the number of particle-pairs (i.e. Cooper pairs) for the state 
1$), i.e. 


N\<^>) = M\<^>). (7) 

The exact solution was reported in [2], and subsequently shown that it could be derived through a 
variety of means including use of the classical Yang-Baxter equation [3], the Quantum Inverse Scattering 
Method [4], or the Gaudin algebra [5]. A second form of exact solution can be obtained by a particle-hole 
transformation, denoted T, which can be defined by 


TNii: = I-Ni, 

T6,T = b], 

T5|T = 6,, 

where T = and I denotes the identity operator. The action of T naturally extends to states. In 
particular if jy) denotes the completely filled state of L particle-pairs then 

lx) = T|0) 


and if (7) holds true then 


NT\<^>) = (L-M)T|$). 


We then find that 

L 

TH{G)T = -H{-G) + 


This provides a simple relationship between the spectrum of the attractive model with G > 0 and that of 
the repulsive model with G < 0. 

Below, we will investigate in more detail the consequences of (8) in terms of the eigenstates. First note 
that it follows that a second exact solution exists whereby for each solution of the coupled equations 


-G- 


2P-L-1 


Vk 


E — 


= E: 


k = i,...,p 


(9) 
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there is an eigenstate of (2) with energy eigenvalue given by 


E = J2zf + iG-l)J2yk, ( 10 ) 

1=1 k=l 

where P = L — M is the quantum number which denotes the number of hole-pairs. In this instance the 
eigenstate is of the form 


where 


3 = 1 


B{y) = TC{y)T 

L 

— _ 5 - 

3 = 1 ^ J 


( 11 ) 


This result can also be obtained by a direct calculation following the methods of [19]. The details are 
provided in Appendix A. 


3 Duality 

Previous studies [4,5] have identified an exact duality relation in the spectrum of the Hamiltonian. Here 
we will first recall the essential aspects of that result, before continuing to establish a second exact duality. 
A direct commutator calculation leads to the result 

[H{G), C(0)] = (1 + G) [Ho, G(0)] - GQl [Q, G(0)] 

= Q\GL-2GN -G-l). 

Using proof by induction we have more generally 

[H{G), (G(0))‘^] = JQ\G{Q)y-^{GL - 2GN -GJ-1) (12) 

At this point it is important to make the observation that TV is a conserved operator, which partitions the 
space of states according to it eigenvalues M = 0,1, 2,. Consider a state |4'') satisfying 

H{G)m = E{G)m, 

= M'l'S'). 

Choosing G~^ = L — 2M — J and setting 

[4/) = (G(0))^|4'') 

it follows from (12) that 

H{Gm = E{Gm, 
iV|4') = (M'-h J)|4'). 

We call jdt) and [T') dual states. Setting M = M' + J dual states are characterised by the relation 

M + M' = L-G-\ (13) 

We may view the state jdt) as a quasi-condensate containing J particle-pairs, each of zero energy in accor¬ 
dance with the expression (6). 
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In the hole-pair picture the duality relation (13) is maintained, and can be verified in a similar fashion. 
In this instance it is convenient to use (4) to express (2) as 

L 

// = G ^ J + (1 - G)Ho - GQQ^. 


It is then found that 


[H{G), i?(0)] = (1 - G) [i?o, i?(0)] - GQ [Qt, ij(o)] 

= Q(2GIV-GL-G-f 1). 

Using proof by induction we have, more generally, 

[H{G), iB{0)y] = JQ{B{0)y-^{2GN -GL-GJ+ 1). (14) 

Consider a state |$') satisfying 

H{G)m = E{G)W), 

IV|$') = 

Choosing G~^ = L — 2M + J and setting 

1$) = (5(0))^!$') (15) 

it follows from (14) that 

H{Gm = E{Gm, 

fV|$) = (M' - J)|4>). 

Setting M = M' — J the states |$) and |<I)') are characterised by the same duality relation (13). Here we 
may view the state |<I>) as a quasi-condensate containing J hole-pairs, each of zero energy in accordance 
with the expression (10). 

3.1 A mixed duality 

Next we turn the discussion toward a different form of duality relating eigenstates of the Hamiltonian H{G) 
with eigenstates of H{—G). We first note the preliminary lemma 

[Ho,Q^] = KQ^-^[Ho,Q] (16) 

which is proved by induction. The calculation is straightforward so we omit the details. Expanding out the 
commutators in (16), this expression can be rewritten as 

HoQ^ - Q^Ho = KQ^-^HoQ - KQ^Hq 
HoQ^ + {K- l)Q^Ho = KQ^-^HoQ. (17) 

Next we proceed to the following identity, which is also proved by induction: 

L 

[Ho, Q^] - [Qt Q^]Q = -K{HoQ^ + Q^Ho) + (18) 

1=1 


The case K = 1 may be verified directly using (4): 


L 

[Ho, Q] — [QI, Q]Q = —{HqQ + QHo) + ^“jQ- 

j=i 


(19) 
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Assuming that (18) holds true for AT = fc — 1, and using (19), we obtain 
[Ho, Q"] - [Q^ Q’^]Q = [Ho, Q’^-^]Q + Q'^-^Wo, Q] 

-{k - l){HoQ'^-^ + Q’^-^Ho) + {k-l)J2 j Q 

+ Q"-' |^-(i?oQ + QHo) + J2 

L 

= (1 - k)HoQ'^ - kQ^-^HoQ - Q^Ho + 

1=1 

Appealing to (17) with K = k then yields 

L 

[Ho, Q^] - [Q\ Q"]Q = -k{HoQ^ + Q^Ho) + kJ2^jQ 

i=i 

which establishes that (18) holds true for K = k. The result follows by induction. 

Finally, 

H{G)Q^ + Q^Hi-G) = ((1 + G)Ho - GQ^Q) ((1 - G)Ho + GQ^Q) 

= HoQ^ + Q^Ho + G[Ho, Q^] - G[Q^Q, Q^] 

= HoQ^ + Q^Ho + G[Ho, Q^] - G[Q\ Q^]Q. 

Through use of (18), this leads to the following identity: 

L 

H{G)Q^ + Q^H{-G) = (1 - KG){HoQ^ + Ho) + KG^, z]Q^. 

1=1 



Consider a state |0') satisfying 

i7(-G)|0') = A'|0'), 

Ar|0') = M'\e'). 

Choosing G~^ = K and setting 

|0) = g^|0') 

it follows from (20) that 

i7(G)|0)=A|0), 
Af|0) = M|0). 

where 

E = j^z^-E', 
1=1 

M = M' - K. 


We call |0) and |0') mixed dual states, which are characterised by the relation 

M' -M = G-^. 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 
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Note that 


Q = lim yB{y). 

y—¥oo 

In contrast to the earlier discussed duality in terms of zero energy hole-pairs, the operator Q is associated 
with an infinite energy hole-pair. However the energy of the state remains finite, as the sum of the diverging 
energies is finite. We prove this in Sect. 4 through an analysis of the Bethe Ansatz equations. 

3.2 The combined duality 

The two dualities described above can be combined to give the following result. Consider a state |H') 
satisfying 


iJ(-G)IH') = A'lH), 
N\n') = M'\n'). 


Setting 


|H) = Q^(i?(0))^|H') 


(24) 


where 


K = G-\ 

J = 2M' - - L, 


then 

H{G)\n) = E\n), 
N\n) = M|0), 

where 


E = 




-E', 


M = M' - K - J. 


The last relation is equivalent to 


M + M' = L. 


(25) 


One feature of the combined duality is that the dimension d{M) of the sector for fixed M is equal to 
d{M'). Explicitly, d{M) is given by the binomial coefficient 


d{M) 


LI 


(26) 


It follows from (25) that d{M) = d{M'). This points to the possibility that the mapping from states |H') 
to |0) through (24) is a bijection, which we believe to be true. Numerical results which support this view 
are discussed in Sect. 5. 


4 Compatibility of dualities and Bethe Ansatz equations 

In [4] it was shown that the duality characterised by (13) is compatible with the Bethe Ansatz solution (5). 
Here we extend that analysis to accommodate the second Bethe Ansatz solution (9). 
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Consider a generic splitting of the set of roots Y of (9), with |y| = P, into non-intersecting sets Y' and 
Z such that Y = F' U Z. We can express (9) as 


-G-^+2P-L-i+y y 

7/7 — r. • ^ 


‘^Uk 


E 


2j/fc 


y^zy^-y^ 


-G-^+2P-L-1+Y' „ = y 

7/7 — r. • ^ 


‘2‘yk 


E 


2yfc 


y.wyj^-y^ ytTv'y^-y^ 

Setting \Y'\ = S and |Z| = T we take the sum in (28) over elements in Z to give 

L 

Vk 


r(-G-‘+2P-L-l)+ E 


‘^Uk 


E 


2/fe Gi"' 
yk G Z. 

2j/fe 


(27) 

(28) 


i/feGZ/=l^'= yj,yk(iZ,yj^yky^ yjeY',ykeZ 


= r(T-i)+ y 


Vj&Y',yk&z 


2yk 

yk - yj' 


(29) 


Suppose that at some limiting value of G we have 

yk^O for all yk G Y, 
yk = 0 for all yk G Z. 

Taking note that S + T = P Eq. (29) informs us that 

T= -G-^ + 2P - L, 

while from (27) we obtain 

L 

-G-^+2P-L-1 + J2 ^’" 2 = E 


(30) 


‘^Ur) 


k^l 

L 


-G~^-b 2 (p-r)-L-i + y = y 


2T, 


ym G F , 


G F'. (31) 


Eq. (31) is the set of Bethe Ansatz equations ior S = P — T roots of an attractive system. The above 
calculations indicate that given such a solution set Y', we can augment it with T additional roots which all 
have zero value to obtain the solution set F, provided that T is given by (30). Identifying 

M = L-P, 

M' = L-S, 

then (30) is equivalent to the duality relation (13). 

Alternatively, suppose that at some limiting value of G we have 


Eq. (29) in this instance leads to 


while from (27) we obtain 


y^^oo ior all ym€Y', 
ym = 00 for all ym G Z. 

T = -G"^ -t 2P - 2S, 


(32) 


-G~^-f2P-L-1-f-y y"'o = y 


‘^Ur) 


k=l -- "fc y.GY'M^yr. y^^ y^ 


G-i+2(p-G-i)-L-i+y-^= y 


2^77 


- Vm , Vm Vj 


y-m ^ E , 
Vm e r. 


(33) 



Eq. (33) is the set of Bethe Ansatz equations for S' = P — G~^ roots of a repulsive system. The above 
calculations indicate that given such a solution set T', we can augment it with T additional roots which 
all have infinite value to obtain the solution set Y, provided that T is given by (32), which simplifies to 
T = G~^. Identifying as before 


M = L-P, 
M' = L-S, 


then (32) is equivalent to the mixed duality relation (23). 
Moreover, from (28) we hnd 


(_G-i + 2P-L-l)y, + ^^^= ^ 


(_G-i+2P-L-1) ^ ^ + 


yk^z vk&z 1=1 


Vk - z[ 


E 


^yl 


E 


, Uk & z 
“^VkiVk - Vj + Vj) 


yk,yj&Z,y^^yk yk&Z,yjGY' 


( G"^+2P 1) E Vk + EE 

Hu — r? ^ 1) E + E 


vkdz vkdz 1=1 


yk&Z yk&Z,yj£Y 


‘^VkVj 
, Vk - Vj 


(34) 


Considering the limit G ^ > T, such that yk ^ oo ior yk € Z, we let 

CT = lim yk- 

O-i^T 

Vk^ Z 

Eq. (34) then reduces to 

L 

J2zl = {l-T-^)a + 2 yj 

k—1 yj^Y' 

and consequently the energy is found from (10): 

L p 

E = Y,zf + ^(G -l)J2yk 

1^1 k^l 

L 

= J2zf + (T-i -1)5] yj + (T-i - l)a 

1=1 yj<^Y' 

= (l+T-i) 5] y,. 

yjdY' 

The above expression is simply (6), with the set {yj} a solution of (5). Proceeding further we have 

E = Y,zf-E' 

Z=1 


in agreement with (22), where 


p' = 5]zf-(i + r-i) ^ 

y,-GV' 


1=1 


is the energy of the dual state for the repulsive Hamiltonian H{—T ^) as given by (10). 
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Figure 1: Distribution of roots for L = 8,P = 3, and ej = j/10. The panel (a) corresponds to the system 
with G = 1/2 + 0.0000001. The inset in panel (a) shows the distribution of roots around Ke{y) = 1. As 
G —>■ 1/2 sixteen roots collapse to the point Re(j/) = 1, as illustrated in panel (b) and its inset for which 
C?= 1/2. 


5 Numerical results 


It has been established that the Bethe Ansatz solution (9) has the property that for particular values 
of G there are roots which are zero, and others which are infinite. To numerically solve (9) near these 
values necessarily means that the elements of the solution set will vary across several orders of magnitude, 
potentially imposing a large computational cost. To alleviate this issue we perform a change of variables: 


Vj 


Zl 


1 + 

1 - n/ 

1 + 

1-e/’ 


such that 


Vj = -1, 

Vj = 1 . 

Under this change of variables the Bethe Ansatz equations (9) become 


Vj = 0 
Vj = oo 


-2G- 


AP-2L-2 L-2P + 2 


1 — v‘^ 

J. ‘^r> 


I-v„ 


+ 


L 

E- 

1^1 


Vn - £l 


P 

E 

q^n 


Vn - Vo 


(35) 

(36) 


(37) 


To numerically solve the above equations we adapt a technique described in [17]. 

The first case we consider is a system with L = 8 and P = 3. In this sector the dimension of the state 
space is 56. Consider G = 1/2. With respect to Eq. (23), the dual sector corresponds to P = 1, which has 
dimension 8. In this instance the mapping between sectors given by (21) cannot be a bijection. 

We perturb the coupling by a small amount and numerically solve (37) for all roots, with the results 
displayed in Fig. 2. It can seen that that there is a subset of roots close to the value 1. As G —1/2 they 
converge to 1. In Table 1 the root sets are sorted according to increasing energy which is computed through 
(10). It is apparent that the eight lowest energy states each have two roots close to the value 1. Between the 
particular sectors P = 3 and P = 1, the operator is an injection in the limit G —>■ 1/2. However since 
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L = 8,P = 3,G= (1/2)-H 0.0000001, and Sj =i/10 

Energy 

Bethe roots 

Energy 

Bethe roots 

-7.54299 

0.99997,1.00003,1.49644 

22.36145 

0.35540,0.54943,0.74665 

1.94946 

0.13030,0.99994,1.00006 

22.43743 

0.74865,0.45717 ±0.04099z 

2.44011 

0.23860,0.99994,1.00006 

22.55229 

0.24377,0.55370,0.74773 

3.08380 

0.34552,0.99994,1.00006 

22.70038 

0.13294,0.55556,0.74847 

3.97611 

0.45217,0.99994,1.00006 

22.92889 

0.24797,0.44046,0.75208 

5.30314 

0.55903,0.99993,1.00007 

23.05279 

0.35360,0.63967 ±0.0410H 

7.49744 

0.66657,0.99992,1.00008 

23.07422 

0.75333,0.32407 ± 0.0348H 

11.87227 

0.77566,0.99989,1.00011 

23.07834 

0.13439,0.44448,0.75282 

12.52145 

0.18033,0.21200,0.93364 

23.12185 

0.24304,0.65007 ± 0.04141z 

12.64936 

0.13716,0.33338,0.93230 

23.20321 

0.13263,0.65541 ±0.04094z 

12.83630 

0.93013,0.31915 ±0.03170z 

23.33978 

0.13784,0.33456,0.75472 

12.87624 

0.13380,0.44202,0.92986 

23.34662 

0.48232,0.53200,0.62467 

12.98062 

0.24692,0.43819,0.92852 

23.49317 

0.18649,0.20946,0.75562 

13.24971 

0.13239,0.55018,0.92567 

23.78106 

0.54736,0.54021 ±0.11610z 

13.31567 

0.92413,0.44485 ± 0.04048z 

24.17544 

0.24905,0.44448,0.63927 

13.34510 

0.24279,0.54830,0.92424 

24.30757 

0.64194,0.33081 ±0.0382H 

13.47474 

0.35396,0.54410,0.92226 

24.32987 

0.13496,0.44873,0.64022 

13.99730 

0.13162,0.65940,0.91663 

24.60177 

0.13850,0.33606,0.64421 

14.07798 

0.24104,0.65843,0.91494 

24.63370 

0.24842,0.51582 ±0.02981z 

14.16919 

0.91099,0.57513 ±0.03533z 

24.73161 

0.13481,0.52553 ±0.03537z 

14.18711 

0.34973,0.65670,0.91260 

24.75667 

0.19766,0.20219,0.64576 

14.34237 

0.45985,0.65262,0.90915 

24.87852 

0.53021,0.35199 ±0.04392z 

17.62521 

0.13125,0.77844,0.85271 

25.18880 

0.42933,0.37736 ±0.09555z 

17.77665 

0.24030,0.77968,0.84490 

25.25408 

0.13948,0.33963,0.53465 

18.00125 

0.34831,0.78269,0.83171 

25.40724 

0.53757,0.20335 ±0.01472z 

18.38397 

0.45676,0.79787 ±0.01289z 

25.57900 

0.13973,0.40092 ±0.00763z 

19.32412 

0.56759,0.77091 ±0.04229z 

25.79494 

0.42898,0.21141 ±0.02721z 

20.41664 

0.67105,0.73644 ±0.10065z 

26.03673 

0.31998,0.22656 ±0.05659z 


Table 1: Numerical solution of the Bethe Ansatz Eqs. (37) with L = 8, P = 3,G = (1/2) + 0.0000001 and 
Ej = j/lO. The energies are calculated through (10) after transforming back to the variables yj and zf 
through (35,36). 


the dimension of the co-domain is larger than the dimension of the domain, not all states in the co-domain 
are of the form (24). 

The second case we consider is a system with L = 8 and P = 6. In this sector the dimension of the 
state space is 28. Again consider G = 1/2. With respect to Eq. (25), the dual sector corresponds to P = 2 
which also has dimension 28. In this instance the mapping between sectors given by (21) may possibly be 
a bijection. 

We again perturb the coupling by a small amount to make it apparent that the solutions we obtain are 
not spurious, and that we can account for the full dimension of the sector. We numerically solve (37), with 
the results displayed in Table 2. The root sets are sorted according to increasing energy which is computed 
through (10). It is apparent that all energy states each have a complex conjugate pair of roots close to the 
value —1, and two real roots close to the value 1. As G 1/2 they collapse to —1 and 1 respectively (results 
not shown). Between these particular sectors P = 6 and P = 2 the operator Q^(B(0))^ is a bijection. All 
states in the co-domain are of the form (24). Although we do not have a proof, we expect that this is true 
in general for sectors related through the combined duality (25). A necessary condition for this to be the 
case is that the mapping (15) is a surjection. This requires that d{M') > d{M) where d{M) denotes the 
dimension formula (26) and M, M' are related through (23). A proof of this result is given in Appendix B. 
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L = 8,P = 6,G = (1/2) -h 0.0000001 and sj = j/10 

Energy 

Bethe roots 

-11.83981209 

-0.99999981 ± 0.00024393/ 0.99997107,1.00002897,1.35407523,9.17378076 

-3.65796172 

-0.99999946 ± 0.00057037/ 0.99996003,1.00004003,1.72872448,0.13265362 

-2.80008720 

-0.99999953 ± 0.00050626/ 0.99995840,1.00004166,1.79893201,0.24145319 

-1.74992090 

-0.99999959 ± 0.00045582/ 0.99995624,1.00004382,1.89428154,0.34849627 

-0.41519412 

-0.99999963 ± 0.00041512/ 0.99995320,1.00004686,2.02752239,0.45498685 

1.37437992 

-0.99999967 ± 0.00038118/ 0.99994857,1.00005151,2.21634859,0.56145363 

4.02553678 

-0.99999970 ± 0.00035212/ 0.99994039,1.00005969,2.48406296,0.66839900 

4.55935661 

-0.99999662 ± 0.00199582/ 0.99993382,1.00006628,0.20680832 ± 0.02043937z 

5.01704811 

-0.99999709 ± 0.00178555/ 0.99993146,1.00006865,0.14110366,0.33689835 

5.89937520 

-0.99999774 ± 0.00148867/ 0.99992647,1.00007365,0.13730723,0.44668753 

6.01657120 

-0.99999807 ± 0.00128240/ 0.99992550,1.00007462,0.33690194 ± 0.03941953z 

6.38806868 

-0.99999819 ± 0.00123967/ 0.99992318,1.00007694,0.25126848,0.44251302 

7.21930929 

-0.99999827 ± 0.00124763/ 0.99991780,1.00008233,0.13570667,0.55544129 

7.70384297 

-0.99999861 ± 0.00104596/ 0.99991363,1.00008651,0.24697377,0.55356595 

8.18565261 

-0.99999896 ± 0.00084940/ 0.99990850,1.00009165,0.46637678 ± 0.03850363z 

8.33050601 

-0.99999890 ± 0.00088637/ 0.99990761,1.00009255,0.35820837,0.54928501 

8.82488231 

-0.99999973 ± 0.00032685/ 0.99992267,1.00007744,2.86124262,0.77673107 

9.40566610 

-0.99999862 ± 0.00107506/ 0.99989980,1.00010037,0.13481277,0.66436225 

9.88812763 

-0.99999887 ± 0.00091253z, 0.99989340,1.00010678,0.24509789,0.66345460 

10.51623429 

-0.99999908 ± 0.00078718/ 0.99988379,1.00011642,0.35409601,0.66184154 

11.36752125 

-0.99999923 ± 0.00069725/ 0.99986777,1.00013248,0.46378746,0.65803096 

11.72101631 

-0.99999930 ± 0.00065253/ 0.99985675,1.00014354,0.59339982 ± 0.01996575z 

13.77136379 

-0.99999885 ± 0.00095049/ 0.99985955,1.00014067,0.13423675,0.77454140 

14.25217295 

-0.99999905 ± 0.00081659/ 0.99984872,1.00015152,0.24402831,0.77416049 

14.88038709 

-0.99999921 ± 0.00071424/ 0.99983265,1.00016762,0.35222296,0.77357060 

15.74513536 

-0.99999932 ± 0.00063925/ 0.99980693,1.00019339,0.46018824,0.77251429 

17.00197794 

-0.99999940 ± 0.00058367z, 0.99976580,1.00023458,0.56880755,0.77001849 

18.42439918 

-0.99999946 ± 0.00054548z, 0.99971075,1.00028967,0.68228619,0.74962317 


Table 2: Numerical solution of the Bethe Ansatz Eqs. (37) with L = 8, P = 6,G = (1/2) + 0.0000001 and 
Ej = j/lO. The energies are calculated through (10) after transforming back to the variables yj and zf 
through (35,36). 


6 Phase diagram 

We introduce the rescaled coupling parameter g = GL and the Hlling fraction x = M/L. It is convenient 
to represent the duality relations in terms of the phase diagram shown in Fig. 2, which depicts six regions 
in the g~^ — x plane. For the attractive model with 5 > 0 the three regions denoted IV, V, and VI have 
been previously identified as the strong pairing, weak pairing, and weak coupling regimes respectively [4]. 
The weak pairing and strong pairing regimes are dual with respect to (13). The boundary between these 
regions is known as the Read-Green line and is given by the relation 

x = \{l-g-^). 

This line extends into the repulsive region g < 0 and provides the boundary between regions II and III 
which are also dual with respect to (13). The boundary between the regions IV and VI is known as the 
Moore-Read line, and is given by 

x = l-g-\ 

The ground state on the Moore-Read line is dual to the vacuum through (13) and has zero energy. 
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With respect to the mixed duality relation (23) regions II and V are dual, as are regions III and IV. 
With respect to the combined duality governed by (25), regions III and V are dual. For a state in the strong 
pairing regime (region V) which may be expressed in the form (24), we define the fraction of zero energy 
hole-pairs as 


ho = 


J 

L 


and the fraction of infinite energy hole-pairs as 


hoo 


K 

T 


where as before P = L — M. The states (24) only exist for certain integer values of G~^, however in the 
thermodynamic limit the values of ho and hoc become dense. The thermodynamic limit is obtained by 
taking the limits 


M —)■ oo, 

L —)■ cx), 

G^O 

such that X and g are finite [4,5]. In this limit we have 

ho = l- 2x- g~^, 

hoo — g 


such that ho + hoc is independent of g. 

If we conduct an analogous analysis in the particle-pair picture, through use of the Bethe Ansatz solution 
(5), we do not obtain a complementary portrait. There are no zero energy particle-pairs, nor are there any 
infinite energy particle-pairs, in the strong pairing regime. This is a key result of this study, that there is a 
clear asymmetry between the hole-pair picture and the particle-pair picture. 


6.1 Inversion 


It is worth briefly mentioning that besides the duality relations discussed above, there exists another type 
of relation which we call inversion. Consider the Bethe Ansatz Eqs. (9) and set Uk = Then 


(-G-1 +2P-L- l)uk + y 


= E 




% - Uk 


_G-i+2P-L-l + y + 


-2 


= E 


2{uj -Uk + Uk) 


= 1 h -Uk Uj-Uk 


-G-i -f 2P -1 + y — = 2P - 2 -f y 


2uk 


Ti^i -Uk 




G-i - 1 

Uk 


tiuk-zi^ .Vfc Uk - Uj 


-G-^ +2P - L-l A 1 .y 

—+ E:—^ = E 


Uk 


1=1 Uk - Zi Uk Uj 


where G ^ = —G ^ +2P — L. If the momentum parameters are chosen such that {zi : I = I, = {z^ ^ : 
I = 1,..., L} inversion maps roots for a Hamiltonian H{G) to a set of roots for the Hamiltonian H{G). It 
provides a invertible mapping between solutions sets in regions I and VI, between solutions sets in regions 
H and IV, while regions HI and V are each stable under inversion. 
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Figure 2: Phase diagram. Regions II and III are dual with respect to the relation (13), as are regions IV 
and V. Regions III and IV are dual with respect to the relation (23) as are regions II and V. These dualities 
combine to give a duality between regions III and V with respect to the relation (25). Our analysis suggests 
that the mappings from II to V and from III to IV through (21) are injections, mappings from III to II 
and IV to V through (15) are surjections, such the the composed mappings from III to V through (24) are 
bijections. 

7 Mean-field approximation 

Our hnal point of discussion concerns a mean-field approximation analysis, which is a standard technique 
applied to the analysis of pairing Hamiltonians in general. It has been previously shown in [4, 5] that the 
mean-held gap and chemical potential equations and the Bethe Ansatz solution (5) for the ground state in 
the continuum limit are equivalent. Here we investigate the extension of that correspondence to include the 
Bethe Ansatz Eqs. (9). 

Using a mean-held approach, in particular where products of operators A and B are approximated as 

AB « A(B} + (A}B - {A){B} 
the Hamiltonian (2) may be approximated by 

n = Ho-^A*Q-^AQ^ + ^-^iiN-M), (38) 

where A = 2G {Q), A = |A|, N is the particle-pair number operator, M = (TV) and /i is a Lagrange multiplier 
which is introduced since the mean-held approximation does not conserve particle number. Setting 

S{z,) = ^(z|-m)^ + z|A2 

the ground-state energy is found to be 

Emin = i ^(z| - m) - ^ ^ 

^ i=i ^ i=i 
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associated to the mean-field ground state 


L 


l^min) — |0) 

i=i 

L 


= Y[iujbj+Vjl) lx) 

i=i 


( 40 ) 


where 




Through use of the Hellmann-Feynman theorem we may take partial derivatives of (38) and (39) to generate 
the following constraint equations: 



L 


L-2M = Y^ 
j=i 


Zj - M 

£{zj) 


(41) 

(42) 


which are known as the gap and chemical potential equations. It is apparent that (41) cannot admit a 
solution when G < 0. However (5) maps to (9) with the change G —>■ —G and changing the quantum 
number to count hole-pairs instead of particle-pairs, while the mean-field wavefunction (40) can be equally 
expressed in terms of particle creation operators acting on the vacuum or particle annihilation operators 
(i.e. hole creation operators) acting on the completely filled particle state. At first sight it appears there is 
a paradox. 

If on the other hand we calculate that highest energy state of the approximation (38) we find 


Ifmax — 2 2 ^ + 

i=i 


4G 


■/iM 


(43) 


associated to the mean-field highest-energy state 

L 

(44) 

= Wivpj - Ujl) lx) • 

i=i 

where * denotes complex conjugation. Through use of the Hellmann-Feynman theorem we may take partial 
derivatives of (38) and (43) to generate the following constraint equations: 


j 


G 

2M-L=y i, ; . 


(45) 

(46) 


Eqs. (45,46) are equivalent to (41,42) via the particle-hole transformation and the change G —G. This 
basically asserts that the mean-field approach is justihed in calculating the low energy spectrum of the 
attractive model or the high energy spectrum of the repulsive model. The observation is entirely consistent 
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with (8). Note that in either case we can project the mean-field states (40,44) onto the sector with fixed 
M, which leads to the following unnormalised states 

M 


M 


It is known [4] that (47) is exactly the ground state on the boundary between regions IV and VI (Moore- 
Read line). Both forms shown in (47) are obtainable by the Bethe Ansatz solutions (5) and (9) respectively. 
By the same methods it can be shown that (48) is exactly the highest energy state on the boundary between 
regions I and II. These boundary lines provide the most “mean-field-like” states. In the hole-pair picture 
they arise from mappings of the form (21) where the domain is the one-dimensional space with basis { 
chi)}. 

What is not apparent is how a mean-field approach may be implemented to observe the structure of 
states with the form (24) in region V, even at the level of an approximation. Conversely there are well- 
established methods (e.g. see [18]) which in principle permit the calculation of the ground state energy 
of the attractive system from the Bethe Ansatz solution (9) in the continuum limit. Whether or not this 
approach simply reproduces the continuum limit of (41,42), or produces some new insights, presents an 
interesting open question. 


i,i=i 


l^'max) -> 


j=l 



8 Conclusion 

We presented an alternative form of the Bethe Ansatz equations, based on the hole-pair picture, in order to 
re-examine the p + ip model. One of the main results we discovered is an intrinsic asymmetry between the 
particle-pair and the hole-pair perspectives, in contrast to the s-wave paring model on which such a symmetry 
can be imposed [13,14]. In particular for the attractive pairing system there are instances of diverging roots 
of the Bethe Ansatz equations in the hole-pair picture which can be precisely identified and counted. It 
led us to conjecture that all states in the strong pairing regime have the form of a quasi-condensate with 
the same number of zero energy pairs, and infinite energy pairs whose energy sum is finite. Significantly, 
diverging roots do not occur in the Bethe Ansatz solution of the attractive model in the particle-pair picture. 
Our findings are summarised in the phase diagram Fig. 2. A notable feature of the phase diagram is that 
the boundary lines, which were determined by exact calculation without approximation, and all independent 
of the parameters zi implicit in the Hamiltonian (2). 
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Appendix A - Direct calculation of the exact solution 


We start with the observation that 


H\x) ='^z‘^\x)- 

i=i 
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To determine exact eigenstates of the Hamiltonian by way of a Bethe Ansatz solution, we follow the approach 
of [19]. Define generic states of the form 

M M M 

IV') = n ^iyk)\x), li’j) = n ^iyk)\x), \^ji) = n ^iyk)\x) 

k=l k^j k^j,l 

where B{y) is given by (11). Noting the commutation relations 


[i/o, 5(2/)] =-E 


^ 4 

^ b 


2^3 


,=,y-^3 

= Q- yB{y) 

[Q\B{y)]=J2-^m,-I) 
i=i y 

it is then found that 

- E = ((1 - + (G - 1 ) ^ ^|/)|vl') - GQQ^\^) 


1=1 


1=1 


= (1 - G) ^ (Qldi,) - y,m -GQJ2J2 (B{yi)- - I)) ...5(yp)) |x) 

,=i ,=ip=i^ \y h / / 


P 




= (1 - G) ^ (gim,) - y.id/)) + Gg^ ^ 


PPL 


+ 2GgEEE 


1=1 r>j p 


^ ( 2/1 - 2p)(2/r- - Zp) 


2\ bpl'i^rj) 


PL 0 

Zf 


= {i-G)j2 (01^.) - 2/,i'i^)) - GgE E ;7E72 

1=1 1=1 1=1 ' 

^ ^ ^ / y z'^ 

^ ^ \ (2/j - yr)iyr - zD ^ (jjr - 2/j ) (2/j - Z^) 


yjh 




L P 2 

zf 


= (1 - G) ^ (gjd/,) - y, Id/)) - Gg ^ ^ ^^jvi/,) 

3=1 1 = 13=1 ' 

p p 


cgEE 




I'i'r) 


= lt^j\yO-yr yr-y3 

= (1 - G) (gid/,) - y,jd/)) -Gqj^jz + 2GgEE 


1=1 




1=1 -r^j 


yj - yr 


The terms proportional to jd/j) cancel provided 

G-I + GV 2 = 2 Gy y^' , k = l,...,P 

S ^ % - yr 

which can be equivalently written as (9). For each solution of those coupled equations, |^/>) is an eigenstate 
of the Hamiltonian with energy eigenvalue given by 

L p 

E = Y^z^ + {G-l)Y,yu. 

Z =1 
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Appendix B - Proof of an inequality 

Here we show that when M + M' = L — G~^ with G~^ > 0, and L > M' > M, that 

d{M') > d{M) 


where d{M) is given by (26). Noting that 

L>M' >L-M'-G-^, 

this is equivalent, in view of (26), to establishing that ior A — G > B >G, 

C\{A-G)\ ^ 

B\{A-B)\ - ■ 


When X > F we have from the definition of the factorial that 

X! 


Then 


This establishes that 


F! 

X! 

F! 


< X 


X-Y 


> (F + 1) 


X-Y 


- c)\ 

B\ 

G\ 


>{B + l) 


A-B-C 


> 


iA-B)\ - {A-B)^-b-c- 
C\{A-C)\ 


B\{A-B)\ 


> 1 


whenever 2B + 1 > A. Alternatively, 

(^ - C)l 

{A-B)\ 

C! 

'W\ 

This establishes that 

C\{A-G)\ 
B\{A-B)\ - 

whenever A + 1 > 2B. It then follows that (49) is true. 


> (A-H + 1)^-^, 
1 

- 


( 49 ) 
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